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SUMMARY 

Results are presented of a theoretical and experimental investiga- 
tion of heat transfer involving laminar natural convection of fluids 
enclosed between parallel walls oriented in the direction of the body 
force, where one wall is heated uniformly and the other is cooled uni- 
formly. For the experimental work, parallel walls were simulated by 
using an annulus with an inner-to- outer diameter ratio near 1. 

The results of the theoretical investigation are presented in the 
form of equations for the velocity and temperature profiles and the ratio 
of actual temperature drop across the fluid to the temperature drop for 
pure conduction. No experimental measurements were made of the velocity 
and temperature profiles, but the experimental results are compared with 
theory on the basis of the ratio of the actual temperature drop to the 
temperature drop for pure conduction. Good agreement was obtained be- 
tween theory and experiment for axial temperature gradients of 40 F per 
foot or larger. 


INTRODUCTION 

Increased application of heat transfer to and from fluids in chan- 
nels has recently required further knowledge as to the heat-transfer 
coefficients and temperature profiles occurring with natural convection. 
Turbine-blade and nuclear-reactor cooling are two of the fields concerned 
with this problem. Work on free-convection heat transfer over a vertical 
plate gave good agreement between theory and experiment. Few results 
have been obtained for the similar case of flow in channels. Reference 1 
obtains a theoretical solution for free-convection heat transfer for 
fluids enclosed in channels . The reference uses a postulated velocity 
distribution to obtain a solution. References 2 and 3 extend this ana- 
lytical work to give an exact solution of the equations in more general 
form for constant wall temperature and constant heat flux, respectively. 
Reference 3 includes the effect of forced as well as natural convection. 



2 


MCA TN 3328 


These three references treat the case of infinite channels with the chan- 
nel axis oriented in the direction of the hody force and are subject to 
the same assumptions; namely , two-dimensional laminar flow, uniform 
axial temperature gradient, and constant fluid properties, except that 
the density is allowed to vary in the bouyancy term. 

The purpose of this report is to compare the results of theory and 
experiment as a check on the assumptions involved in the analytical work. 
For mathematical simplicity, the flow between two infinite parallel 
plates was considered with one plate heated uniformly and the other 
cooled uniformly. An exact solution was obtained with the assumptions 
cited previously for references 1 to 3. In order to simulate the case 
of infinite parallel plates experimentally, an annulus formed by two 
concentric tubes was used. For all practical purposes, the walls of 
these tubes can be considered parallel if the ratio of the radii of the 
tubes is near 1. This radius ratio limits the spacing between the con- 
fining walls (and, consequently, the range of Grashof numbers obtainable), 
unless very large diameters are used. 

A comparison between theory and experiment is given in this report 
on a heat-transfer basis alone, since no temperature or velocity profile 
measurements were made. The work was done at the MCA Lewis laboratory. 


AMLYSIS 

Steady-state heat transfer through a fluid enclosed by two infinite 
parallel plates oriented in the direction of the body force is considered. 
One plate is heated uniformly and the other is cooled uniformly. The 
flow is laminar and parallel to the body force. The fluid properties 
are assumed constant, except that the density is allowed to vary in the 
bouyancy term. Viscous dissipation and work against the force field are 
neglected. It is further assumed that the axial temperature gradient is 
a constant throughout the system for any particular set of conditions. 

When the aforementioned conditions are applied, the energy equation 
reduces to 

k S 2 t _ St 
P°p Sx 2 = ^ ^ 

(Symbols are defined in appendix A, and a detailed discussion of 
the analysis is given in appendix B.) 

Similarly, the Navier-Stokes 1 equations reduce to one equation in 
the z -direct ion: 


(Bl) 



(B7) 



MCA TN 3328 


3 


The reference temperature t is taken at the center of the channel 


where the net viscous force is zero 


/ a 2 w \ 

U 2 " 0 j ■ 


For a closed system of unit width with zero net mass through-flow 
and constant density. 


i/T w 


dx = 0 


where x is taken to he zero at the center of the channel. 

From these three equations and the appropriate boundary conditions, 


(1) 

\ / x _ r 


( 


Wx=-r k 


w x=r “ 0 


w x=-r “ 0 


the equations for the velocity and temperature profiles in dimensionless 
form are, respectively. 


W* =s 4v 


. vx , vx 
sin — cosh — 
r r 


. , vx vx 
sinh — cos — 
r r 


sin v cosh v sinh v , cos v 
+ — : + 


sinh v cos v 


sin v 


cosh v, 


(B17) 


VX vx vx vx"' 

sinh v cos v + sinh — cos — cosh v sin v + cosh — sin — 

= _ _±_ | £ £ + £ r 

2v l cosh v . sin v sinh v , cos v 


cos v 


sinh v 


sin v 


cosh v 


(B19) 


where, by definition, 


v 


■V 


'PrGr 

64 


Pr 


V 

k 

2 n St 


Gr = 


g p P s 
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The temperature drop across the fluid from wall to wall is then 
found from equation (B19) when x = r: 


e 


kv 


o 2 

sin v + slnh v 


; ) 


(B20) 


sinh v cosh v + sin v cos v 

i J 

The temperature difference across the fluid for pure conduction is 

( 0 r) c = " (B2l) 


The 
tion is. 


ratio of the actual temperature drop to that for pure conduc- 
therefore. 


CD =5 

^r 


0 r 


(= 


2 2 
sm v + sinh v 


sinh v cosh v + sin v cos vj 


(B22) 


These equations apply equally well to force fields other than gravi- 
tational, if the constant g is replaced by the constant of the force 
field of interest. 


EXPERIMENTAL APPARATUS 

A schematic diagram of the experimental equipment is shown in fig- 
ure 1. In order to simulate infinite parallel plates, an annulus formed 
by two concentric stainless -steel tubes was used to contain the fluid. 

1 3 

The resulting annulus was 10^ inches long with an outside diameter of 1— 

Ci 8 

inches and an inside diameter of 1 inch. Thus, the spacing s between 
the walls was 3/l6 inch. The tube walls were l/l6 inch thick. The outer 
tube was heated with a Nichrome -wire element spiral -wound around the 
tube, with a wire spacing of 3/32 inch, which probably gives fairly uni- 
form heating of the hot wall. The inner tube was cooled by forced air. 
The cooling-air-passage length-to-diameter ratio was 84. A wire in the 
form of a helical spring was inserted in the cooling passage to increase 
the heat-transfer coefficient. This system should provide an essentially 
constant heat -transfer coefficient, which is necessary to provide uniform 
cooling of the cold wall with a linear variation of wall temperature. 

The cooling-air-flow rate was measured with a calibrated rotameter. 

The- electrical power was supplied by a transformer with a variable 
secondary voltage and rated at 1 kva at maximum voltage. The electrical 
power input was measured with a calibrated ammeter and voltmeter. All 
temperatures were measured with iron-constantan thermocouples and a 
self-balancing potentiometer. The thermocouples used to measure wall 
temperatures were set into slots milled axially into the tube walls. 
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The thermocouple leads were enclosed in Inconel tubing laid in t h e slots. 
The thermocouple junctions were covered with insulating cement, and the 
slots were then filled with silver solder and finished off flush with 
the tube surface. 


RESULTS AND DISCUSSION 
Analytical 

The dimensionless velocity and temperature, which can be calculated 
from equations (B17) and (B19), respectively, are shown graphically in 
figures 2 and 3 for a few values of the product of Prandtl times modi- 
fied Grashof numbers . The equations show that the dimensionless velocity 
and temperature depend only on this product and the position in the chan- 
nel. The Grashof number resulting from the analysis differs from the 
conventional Grashof number insofar as the temperature difference usually 
appearing is replaced by the temperature gradient and the characteristic 
length appears as the fourth power instead of the cube. The product of 
Prandtl number times Grashof number is sometimes referred to as the 
Rayleigh number. 

The velocity profiles of figure 2 show the point of maximum velocity 
moving closer to the wall as the Rayleigh number is increased. It is 
felt that this fact, in combination with the stabilizing influence of 
the wall, would tend to inhibit the occurrence of turbulent flow. 

The contribution of convective flow in the heat-transfer process 
can best be seen from the ratio of the actual temperature drop across 
the channel to the temperature drop for pure conduction. This ratio, 
which is the same as the ratio of the molecular to the apparent conduc- 
tivity, can be obtained from equation (B22) and is plotted in figure 4. 
For large values of PrGr, the actual temperature drop across the fluid 
is less than l/3 the temperature drop for pure conduction. 


Experimental 

A typical plot of measured wall temperatures for large temperature 
gradients is shown in figure 5. With uniform heating of the hot wall 
and constant heat -transfer coefficient on the cold wall, the temperature 
gradients on the hot and cold wall are constant and equal, except for 
the extreme ends of the test section. Data were not obtained to verify 
the assumption of -uniform temperature gradient throughout the fluid. 

For small temperature gradients the wall temperatures were linear, 
but the gradient on the heated wall was larger than that on the cooled 
wall. The assumption of uniform temperature gradient throughout the 
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system is, therefore, invalid for small temperature gradients. With 
water as the test fluid and with 3/l6-inch spacing between the plates, 
the m inimum value for which the assumption of uniform temperature gradi- 
ent is valid is approximately 40° F per foot. Additional data are re- 
quired to determine this limit more accurately. The data presented in 
this report are for temperature gradients of 40° F per foot or greater. 

Because of the end effects on the flow and because of the heat 
losses from the ends of the test section, it was impossible to obtain 
a heat balance between the electrical heat input and the enthalpy rise 
of the coolant. Measurements of the radial temperature drop through the 
insulation surrounding the heater indicated a negligible radial heat 
loss. In that section of the tube where the wall temperature varies 
linearly, there is no net axial conduction in the wall. Hence, the uni- 
form heat flux (as calculated by the electrical heat input) that enters 
the wall is transmitted directly to the fluid. The heat loss is limited 
to the ends of the test section and, therefore, does not affect the re- 
sults of the test. 

The calculated temperature drop through the walls containing the 
water was negligible compared with the temperature drop across the water, 
and, therefore, the measured wall temperatures were taken to be the tem- 
perature of the surface adjoining the water. 

The temperature drop across the fluid for pure conduction At c was 
calculated from the equation for an annulus given in reference 4: 

r o l°Se( r oAi> 

At <= " 4 o k 

where q" was calculated from measurements of the electrical heat 
input. 

The physical properties used in evaluating Prandtl number, Grashof 
number, and the temperature drop for pure conduction were obtained from 
reference 5 for saturated liquid water, except for the coefficient of 
thermal expansion (3. Values of the coefficient of thermal expansion 
of water were taken from reference 6 and are plotted in figure 6. These 
values are mean values for 1° temperature change. Inasmuch as p is a 
function of temperature, different mean values would be obtained for 
larger temperature increments. The properties were all evaluated at the 
bulk temperature of the water for each test. For flat plates, the bulk 
temperature is the center-line temperature halfway between the ends. 

Figure 7 shows the ratio of the actual temperature drop across the 
fluid to the temperature drop for pure conduction plotted against the 
product of Prandtl times Grashof numbers. The data points on the figure 
are the experimental results taken with water as the test fluid, and the 
curve is the result of analysis taken from figure 4. 
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Comparison Between Analysis and Experiment 

Inasmuch as no temperature measurements were made within, the fluid, 
a comparison between theory and experiment can be made only for the over- 
all heat-transfer results. Except for the random scatter of the heat- 
transfer data shown in figure 7, good agreement is obtained between 
theory and experiment. Inasmuch as the same assumptions were used in 
references 2 and 3, the agreement shown here lends support to the more 
general analytical treatment given .in these references. 

The effect of the curvature of the walls on the experimental re- 
sults is unkn own. j but this effect is believed to he minimized, because 
the data are in the form of the temperature ratio having this effect in- 
cluded in both the numerator and the denominator. A theoretical inves- 
tigation for an annulus is necessary to permit a comparison between 
theory and experiment with larger diameter ratios for the annulus. 

In order to determine the upper limit for which the analysis is 
valid, it is necessary to obtain data at higher values of the product 
of Prandtl number times modified Grashof number. 


CONCLUSIONS 

For the range of conditions investigated, the simplifying assump- 
tions used in the analysis for free convection in channels oriented in 
the direction of the body force are reasonable and lead to accurate 
quantitative results. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 8, 1954 
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APPEHDIX A 


SYMBOLS 

The following symbols are used in this report; 

A 1- ,A 2 ’ ’ ’ cons tants 

V/SP P c p St ^+.-1 


B 1> B 2‘ 


D 

F 

z 

Gr 

g 

i 

k 

Pr 

P 


r 

r J 


constants 

constant of integration 

specific heat at constant pressure, Btu/(lb) (°F) 
operator, S/Sx 

body force in z-direction, lb/(hr B )(sq ft) 

2 0 St 4 
SP P ^ s 

modified Grashof number, 5 , dimensionless 

acceleration due to gravity, 4.17 KLO® ft/hr^ 

thermal conductivity, Btu/ (hr) (sq ft)(°F/ft) 

Prandtl number, c^p/k, dimensionless 

pressure, lb/sq. ft 

heat flux at wall, Btu/ (hr ) (sq ft) 

heat flux at outer wall of annulus, Btu/(hr) (sq ft) 

l/2 plate spacing, s/2 

inside radius of annulus, ft 

outside radius of annulus, ft 

plate spacing, ft 
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t 


t 


f 


t 


r 


t 

-r 

u 


v 

w 


x 

z 

P 

e 



n 

p 


temperature, °F 
reference temperature, °F 
temperature at x = r, °F 
temperature at x = -r, °F 
ax/-^, dimensionless 


“/V*. - ]/W 


velocity, ft/hr 
dimensionless velocity, 


wpc p s 


St 


transverse coordinate, ft 
longitudinal coordinate, ft 

o - 1 

coefficient of thermal expansion, F 
t - t_ r , °F 



temperature difference across fluid for pure conduction 

(t - t_ r )k 

dimensionless temperature, n 

Q. s 

viscosity, lb/ (hr ) (ft ) 
density, lb/cu ft 

ratio of actual temperature difference across fluid to 
temperature difference for pure conduction 
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APPENDIX B 


DETAILS OF ANALYSIS 


For the conditions of the problem stated in the ANALYSIS, the gen- 
eral energy equation reduces to 


k 


_Sft 

P°p Sx Z 


« W 


St 


(bi) 


Similarly, the Navier-Stokes ' equations for incompressible flow reduce 
to 


Sp 


= 0 


K. - g 


a w _ 
11 Sx 2 " 


(B2) 


The body force under the action of gravity is 

' N 

Since the pressure is independent of x, then -r 
of x and can be evaluated at any value of x. It is convenient, how- 

O W / \ 

= 0. Then, from equations ^B2j and 


ever, to evaluate -ie- where _ 

^ &x 2 


(B3) 

is independent 


(B3), 


= - P f 


(B4) 


a 2 


w 


where the subscript f denotes a reference condition where = 0. 

Sx^ 

Combining equations (B2), (B3), and (B4) gives 
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From the definition of the coefficient of thermal expansion (3, 

> 1 + p(t - tj,) 

Combining (B5) and (B6) yields 


Pf 

P 


(B6) 


H = (t f - 1) 


dx 




(B7) 


The reference temperature tj is a function of z and varies 

linearly with z as does t } according to the assumption of constant 
temperature gradient. The velocity w is therefore only a function of 
x, and its derivative is written as a total derivative. 


Differentiating equation (B7) with respect to x gives 

d 5 w _ gpp St 


dx- 




(B8) 


Differentiating with respect to x again gives 


-4 
d w 


dx' 


.4 


gpP S^t 
* Sx 2 


(B9) 


Combining equations (Bl) and (B9) yields 


,4 
d w 

dx 4 


SP Pc. 


P St 


pk dz 


w 


(BIO) 


4 SP P c p St 
a = jlk Then 

,4 

d w , 4 ~ 

— 7 + a w = 0 

dx 

The solution of this homogeneous equation can he obtained by using 
the operator D: 

(D 4 + a 4 )w - 0 
(D 2 + ia 2 ) (D 2 - ia 2 )w = 0 
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(D + i -\fL a )(B - i-y/i a)(D + -y/i a)(D - -y^i a)w = 0 



i + 1 

V2 


and 



Therefore, 



log w + — ax = C 

V2 


C 

e 


i - 1 

- V2 


ax 

= w 


w 


C " 
e e 



iax ax 

"Vi Vi 

e v 


(B12) 


Similarly, from equation (Bll), 


iax ax 

V2 V2 

; e v 

(B13) 

iax ax 


; V2 e " V2 

(B14) 

iax ax 


V2 V2 

: e 

(B15) 
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The general solution is equal to the sum of equations (B12) to 
(B15), or 

w = e e U + Aj e -U ^ + e^ u ^Ag e -u + A^ e u ^ 

ax 

where u = Z^2' 

Using the Euler formulas , 

-iu 

e = cos u - 1 sin u 

iu ... 

e = cos u + l sin u 

and the conversion formulas 

e u = cosh u + sinh u 

e u = cosh u - sinh u 


and letting 


B 1 ~ A 1 + A 4 + A 3 + A 2 


B 2 ~ A 1 + A 4 " A 3 " A 2 


B 3 ~ i ^ A 2 A 3 + A 4 A l^ 
B 4 - i(A 4 - A x - A 2 + A 3 ) 


the velocity equation Becomes 

w = B-j_ cos u cosh u + Bg cos u sinh u + B^ sin u cosh u + 

B^ sin u sinh u (B16) 

The constants of equation (B16) can he evaluated from the Boundary 
conditions. There are four Boundary conditions required to determine 
the four constants. The velocity at the Bounding walls must Be zero, 
while the Boundary conditions on the temperature are given By the tem- 
perature gradients. It can Be shown from equation (Bl) that, for the 
case of interest here, the temperature gradients at the two walls must 
Be equal. Integrating equation (Bl) gives 
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For a system closed at both ends, there is no net mass through-flow; 
therefore, with constant density. 



w dx = 0 


Consequently, the transverse temperature gradients at the two walls must 
be equal. The boundary conditions are expressed mathematically as 



w = 0 
x=r 



When these boundary conditions are imposed on equations (B8) and (B16), 
the constants of equation (B16) are found to be 



B 2 = 


4q"v sin v cosh v 

St (sin v cos v + cosh V sinh v) 


pc P s SS 


B„ = 


P c ^s 


4q"v cos v sinh v 

St (sin V cos v + cosh v sinh v7 
; p“ 'Si 


b 4 = o 

Substituting these values for the constants of equation (B16) gives the 
equation for the velocity in dimensionless form: 
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w 


.* _ 


4v 


sin — cosh ~ 
r r 

sin y cosh v 

sinh v cos v 


„ . , vx „„„ vx 
sinh — cos — 


sinh v cos v 


(B17) 


sin v 


cosh 


Differentiating equation (B17) and combining with equation (B7) 
result in an equation for the temperature in dimensio nl ess form: 


tk 

q"s 


t f k 

_t! 

q s 



, vx , vx 
sin — cosh — 
r r 

sinh v ^ cos v 

sin v 


VX 

sinh — cos 
r 


sin v 


cosh v sinh v 



(B18) 


The actual temperature cannot be determined from this equation, because 
the reference temperature t f is also a function of z. 


The difference between the wall, temperature and the temperature at 
any point in the fluid at the same value of z can be determined from 
equation (B18). From the definition of 0*, equation (B18) becomes 


0 * - - 

2v 


VX yx VX VXA 

sinh v cos v + sinh — cos — cosh v sin v + cosh sin —~ 

r _r_ r r 


cosh v sin y 
cos v 


sinh v 


sinh v + cos v 
sin v 


cosh v 

(B19) 


The temperature difference wall-to-wall can be found from equation 
(B19) when x = r: 


q - q"s f sin 2 y + sinh 2 v ^ 

r ky l sinh y cosh v + sin v cos v> 


(B20) 


The temperature difference across the fluid for pure conduction is 


(O « - 


q"s 


(B2l) 


The ratio of the actual temperature drop to the teurperature drop 
for pure conduction is 


9 , 


<Pr " T 


2 2 
sin v + sinh v 


sinh y cosh v + sin v cos vj 


(B22) 
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From the definitions of v and a, 

2p | s4 Vv\ GrPr 

V? y\ k /~ 64 
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Figure 2. - Dimensionless velocity profiles for various values of product of Prandtl number 
times modified Grashof number. 



Dimensionless temperature 



Figure 3. - Dimensionless temperature profiles for various values of product of Prandtl number times 
modified G-rashof number. 
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Prandtl times Grashof number, PrGr 


Figure 4. - Theoretical variation of ratio of actual temperature difference across fluid to temperature difference 
for pure conduction with product of Prandtl times modified Grashof numbers . 




wall temperi 





Temperature, t, °F 


Figure 6. - Variation of coefficient of thermal expansion of water 
with temperature for 1° change in temperature. 
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O Experimental data with water in 



annulus (— = 1.375\ 

\ r i J 



Flat-plate theorv J 





Prandtl times Grashof number, PrGr 

Figure 7 . - Variation of ratio of actual temperature difference across fluid to temperature difference for pure conduction 
with product of Prandtl times modified Grashof numbers . 



